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We investigate non-linear Schrodinger-type formulation of cosmology of which our cosmological 
system is a general relativistic FRLW universe containing canonical scalar field under arbitrary 
potential and a barotropic fluid with arbitrary spatial curvatures. We extend the formulation to 
include phantom field case and we have found that Schrodinger wave function in this formulation 
is generally non-normalizable. Assuming power-law expansion, a ~ t'^ , we obtain scalar field po- 
tential as function of time. The corresponding quantities in Schrodinger-type formulation such as 
, Schrodinger total energy, Schrodinger potential and wave function are also presented. 
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C/^ ; I. INTRODUCTION 
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Canonical scalar field plays important role in inflationary phase in the early universe as well as acceleration in 
' q'' the late universe observed and confirmed by cosmic microwave background large scale structure surveys and 
supernovae type la Q. The scalar field is considered as inflaton field in inflationary models [J]. It could also be 
^ r considered as dark energy that drives the late acceleration described in review literatures 5] and references therein. 

In standard cosmology with Friedmann-Lemaitre-Robertson- Walker (FLRW) background, major components of the 
I— ' late universe are mixture of dark matter which is a type of barotropic fliud and dark energy in form of scalar field. 
When assuming pure scalar fluid in fiat universe, one can obtain analytical solutions otherwise the problem can also 
be solved numerically. However, considering arbitrary types of barotropic fluid and a non-flat universe, it is not always 
' possible to solve the system analytically. 
, Apart from standard cosmological equations, there are few alternative mathematical formulations which arc also 
\^ • equivalent to the scalar field cosmology with barotropic fiuid. One is in form of non-linear Ermakov-attemping equation 
' [g| and another idea proposed recently is in form of non-Ermakov-Milne-Pinney (non-EMP) equation. Cosmological 
QQ , equations in the latter proposal can be written in form of a non-linear Schrodinger-type equation when imposing 
' relations between quantities in standard cosmological equations and Schrodinger-type equation [7] . In case of Bianchi 
, I scalar field cosmology, recent work shows that it is possible to construct a corresponding linear Schrodinger-type 
' equation by redefining cosmological quantities 8] . With the new representation, scalar field cosmology is reinterpreted 
^ ' in new way which might be able to give new methods of approaching mathematical problems in scalar field cosmology. 

There are various observations allowing scalar field equation of state coefficient, to be less than -1 [9|]. Recent 
data such as a combined WMAP, LSS and SN type la without assuming flat universe, puts a strong constraint, 
. W(j, = — l.OBto^os 0- ^Iso the first result from ESSENCE Supernova Survey la combined with SuperNova Legacy 
Survey la assuming flat universe, gives a constraint of — — 1.07± 0.09 [ll|. Therefore it is possible that the scalar 
field dark energy could be phantom, i.e < — 1 [l^ . The phantom behavior, w^j, < — 1 can be attained by negative 
kinetic energy term of the scalar field density and pressure. In FLRW standard cosmology, the field can yield big rip 
singularity, i.e. a, p, \p\ — > cxd at finite time [13| with attempts of singularity avoidance in several ways |14| . 

In this work, we investigate connection between standard cosmological equations and non-linear Schrodinger-type 
equation with a comment on normalization of the wave function. We modify the work of [7] to include phantom field 
case. A case of power-law expansion with scalar field and dark matter is considered as a toy model. We begin from Sec. 
im where we introduce our cosmological system. Afterward in Sec. IIIIl we discuss how non-linear Schrodinger-type 
formulation quantities are related to quantities in standard scalar field cosmology. In non-linear Schrodinger-type 
equation, one important quantity is wave function. We comment on normalization properties of the wave function 
in Sec. IIVI We consider a case of power-law expansion in Sec. |V| before deriving scalar field potential, Schrodinger 
potential and Schrodinger wave function. At last we conclude this work in Sec. IVII 
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II. COSMOLOGICAL EQUATIONS 

In a Friedmann-Lemaitre-Robertson- Walker universe, the Einstein field equations are 
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where = SttG = G is Newton's gravitational constant, Mp is reduced Planck mass, k is spatial curvature, 

Pt and Pt are total density and total pressure, i.e. pt = p^ + p^i and pt = p^ + p^. The barotropic component is 
denoted by 7, while for scalar field, by </>. Equations of state for barotropic fluid and scalar field are p~^ = 



P4> = W(j,P(j,. We consider minimally couple scalar field with Lagrangian density, 

£ = - V{4>) , 

where e = 1 for non-phantom case and —1 for phantom case. Density and pressure of the field are given as 



therefore 



The field obeys conservation equation 
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For the barotropic fluid, we set = (n — 3)/3 so that n = 3(1 + w-y). Hence for cosmological constant n = 0, for 
fluid at acceleration bound {w^ = —1/3) n = 2, for dust n = 3, for radiation n = 4, and for stiff fluid n = 6. Solution 
of conservation equation for a barotropic fluid can be obtained directly by solving the conservation equation. The 
solution is 
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where a proportional constant D > 0. Using Eqs. ([T]), (jl]), ([5]), ^ and ([8]), it is straightforward to show that 
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Therefore if one knows how the scale factor evolves with time, the scalar field velocity and potential can always be 
expressed as a function of time explicitly. 



III. NON-LINEAR SCHRODINGER-TYPE EQUATION 



Non-linear Schrodinger-type equation corresponding to canonical scalar field cosmology with barotropic fluid is 
given by 

d2 



da; 



^u{x) + [E - P{x)]u{x) 
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Quantities in the Schrodinger-type equation above, e.g. wave function u{x), total energy E and Schrodinger potential 
P{x) are related to the standard cosmology quantities as 

u{x) EE a{ty/^, (13) 

2 2 

E ^ -^D, (14) 

P{x) EE !^a{tre^{tr. (15) 

The mapping from cosmic time t to the variable x is via 

x = a{t), (16) 

such that 

ait) = uix), (17) 
m = V'(a;). (18) 



We notice that relation 
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V''(a;)2 = —P{x) (19) 



in Ref. [3] which gives = ±(2/kv^) / V-^i^) '^^ does not include phantom field case. In order to include the 
phantom field case, we modify relation ip{t) = i-0'(x) in Q to e</>(i)^ = x^tip'{xY of which the field kinetic term 
{(fy^) is considered instead of the field velocity (0) so that the parameter e can be included. Therefore, to include the 
phantom field case, corrected relation to Eq. is 



ei^'{xf = ^P{x), (20) 



and i^ix) should read 



K\/n 



V'(x)=±^ [J^dx. (21) 



Inverse function of ipix) exists if P{x) ^ and n ^ 0. It is important for ijj ^(x) to exist as a function since existence 
of the relation x = a{t) (Eq. needs a condition, 

X = Vj"^ o 0(t) = a(t) . (22) 

In case that P{x) = and n 0, then tp = C, hence inverse of ip is not a function since one x gives infinite values of 
-0^^. In this case the relation ((22)) is invalid. If the inverse function, 0;"^ exists (i.e. P{x) ^ and n ^ 0), then the 
scalar field potential, V o a^^(x) can be expressed as a function of time, 

, 12 /du\^ 12^2 3fcu4/n 
^W^^hr -— + + (23) 

Although the potential obtained is not expressed as function of (p, however if one can integrate Eq. PH)) to obtain 0(t), 
the obtained solution can be inserted into a known function V{(j)) motivated from fundamental physics. Then one can 
check which fundamental theories give a matched potential to V{t). The Eqs. pT|) and ([^5)) are indeed equivalent. 
Both require only the knowledge of a(t), D and k which can be constrained by observation. Therefore V{t) in both 
Eqs. ([TT1) and (P^)) can be constructed if knowing these observed parameters. To construct V{t) in Eq. one 
needs to know a{t) as a function of time in order to find u{x) and P{x). However, in constructing V{t) in Eq. (jlip . 
if knowing a(t), D and k, one can directly use these quantities without employing Schrodinger-type quantities. 

IV. NORMALIZATION CONDITION OF WAVE FUNCTION 

Normalization condition for a wave function u(x) in quantum mechanics is 

\u{x)\^dx = l. (24) 
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The wave function here expressed as u[x) = a "/^ — i{t), when applying to the normahzation condition, reads 

/oo 
x^dx = 1 . (25) 
-oo 

In order to satisfy the condition, x must be constant and so is t. Since the form of the wave function must be 
u{x) — i{t) in order to connect equations of cosmology to the Schrodinger-type formulation, therefore u{x) as defined 
is, in general, non-normalizable. 



V. POWER-LAW EXPANSION 



Here in this section, we apply the method above to the power-law expansion in scalar field cosmology with barotropic 
fluid in a non-flat universe. The power-law expansion of the universe during inflation era, 

a{t) = , (26) 

with q> \ was proposed by Lucchin and Matarrese to give exponential potential 

exp|-Ky|[0(t)~0(io)]| , (27) 

assuming domination of scalar field, negligible radiation density and negligible spatial curvature. Recent results from 
X-Ray gas of galaxy clusters put a constraint of g ~ 2.3 for fc = 0, g ~ 1.14 for k = —1 and q ^ 0.95 for fc = 1 [l6j . 
Considering mixture of both fluids, we use effective equation of state, Wes = {p^w^ + p^w-y)/ pt- For a flat universe, 
the power law expansion, a = f^, is attained when —1 < < —1/3 where q — 2/[3(l -I- Wcff)]. If using q ~ 2.3 as 
mentioned above, it gives Woff = —0.71. 



y(</)) = 



g(3g-l) 



A. Relating Schrodinger quantities to standard cosmological quantities 

Assuming power- law expansion and using Eqs. and p7|) . Schrodinger wave function is related to standard 

cosmological quantity as 

u(x) = a{t) = t-'?"/^ . (28) 
We can integrate the equation above so that the Schrodinger scale, x is related to cosmic time scale, t as 

p 

where /3 = {qn — 2)/2 and t is an integrating constant. The parameters x and t have the same dimension since /3 is 
only a number. Using Eq. (pS)) . we can find e4>{t)^ from Eq. ^ 



•/ n9 2q 2fc nZ? 



We use Eqs. (|26p and ([50)1 in Eq. therefore the Schrodinger potential is found to be 



p(,) ^ ^ _ ^ (31) 



With E = —K n D/l2, the Schrodinger kinetic energy is 
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FIG. 1: Potential V{t) plots from non-linear Schrodinger-type formulation assuming a ^ t'', q = 2.3 in flat universe (fc = 0). 
The thickest line is when there is no barotropic fluid D = 0. The middle line is when there is dust fluid together with scalar 
field, i.e. D ^ and n = 3. The small line is when the universe has scalar field with radiation fiuid, i.e. D and n — 4. We 
set K = 1 and in the last two plots, we set D = 1. All plots match results obtained from standard cosmological equations. 



B. Scalar field potential V{t) 
In order to obtain V{t) in Eq. (pS)) . we need to know derivative of u{x): 



= (33) 
At this step, using Eqs. ^ and ^ in Eq. we finally obtain 

Assuming flat universe (fc — 0) and q = 2.3, we show V{t) in Fig. [1] Thickest line on top is of the case scalar field 
without barotropic fluid. The middle line is the case when the dust is presented with scalar field {D ^ 0,n = 3). The 
bottom line is the case of radiation (D ^ 0,n = 4). The V{t) plots from the Schrodinger-type formulation matches 
the plots from standard cosmological equations. The result is independent of e values. The solution (p{t) of Eq. ([50]) 
can not be integrated if e = — 1 or if the integrand of Eq. ((30)) is imaginary. When e = 1 with dust {D ^ 0,n — 3) 
and q = 2.3, the integrand is imaginary. We therefore assume g = 2 to show numerical integrations in Fig. [2] for the 
case D = 0, k = and the case D Q,n = 3, k — 0. In the pure scalar field case D = 0, k = 0, numerical solution 
matches the analytical solution (f>{t) — {y/2q/K)\n{t). This solution can be substituted into Eq. ([M)) to obtain Eq. 
(|27| as in [15j (setting = 1 and 0(io) = 0). When considering cases of closed, flat and open universe containing 
dust matter, V{t) of each case is presented in Fig. [3] where q — 2 is assumed in all cases so that we can see how the 
plots change their shapes when k is varied. It is worth noting that reconstruction of scalar field potential assuming 



scaling solution was considered before in [17 1 



C. Sciirodinger potential P{x) 



We can find Schrodinger potential P{x) from Eqs. (j29|l and ([3T|) where time is expressed as a function of x as 

1 



tix) 



(35) 



Therefore 
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2q(n-2)/(g«-2) ^2^2^ 
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FIG. 2: 0(f) for power-law expansion a ^ f , 5 = 2 in flat universe (fc = 0). The red line is of the when the barotropic fluid 
density is negligible. The green line is in the presence of scalar field with dust (D 7^ and n = 3). In the figure, k = 1 and 
D = l. 



As in Eq. ([5^ . the Schrodinger kinetic energy is 



T(x) 



2qn 



(qn — 2)2 {x 



kn 

T 



{qn — 2) (a; — r) 



2q(n-2)/(gri-2) 



(37) 



The kinetic term has contribution only from the power q and spatial curvature k. A disadvantage of Eq. (|36p is that 
we can not use it in the case of scalar field domination as in inflationary era. Dropping D term in Eq. (j36p can not 
be considered as scalar field domination case since the barotropic fluid coefficient n still appears in the other terms. 
The non-linear Schrodinger-type formulation is therefore suitable when there are both scalar field and a barotropic 
fluid together such as the situation when dark matter and scalar field dark energy live together in the late universe. 
The Schrodinger potentials P{x) plotted with x for power-law expansion with q ~ 2 in closed, flat and open universe 
are shown in Fig. [H In the figure, the dust cases are shown on the right and radiation cases are on the left. We set 
K — 1, D — 1 and r = 0. 



D. Schrodinger wave function u{x) 



The quantity analogous to Schrodinger wave function can be directly found from Eqs. 



and (I55l) as 



u{x) = 



--qn + 1 ] (x-t) 



1 9"/(<3n-2) 



(38) 



which is independent of the spatial curvature k or the initial density D. However, coefficient n of the barotropic fluid 
equation of state and q must be expressed. Wave functions for a range of barotropic fluids are presented in Fig. [S] 
The result is confirmed by substituting Eq. ([55]) into Eq. (H^). 



VI. CONCLUSIONS AND COMMENTS 



We consider Schrodinger-type formulation for a system of canonical scalar field and a barotropic fluid in standard 
FLRW cosmology with zero or non-zero spatial curvature. In the Schrodinger-type formulation, all quantities in 
cosmology are represented in Schrodinger-like quantities and the equation relating these Schrodinger-like quantities is 
written as a non- linear Schrodinger-type equation. If a{t) is known as an exact function of time, a connection of two 
scale quantities, x and t can be found and then other Schrodinger-like quantities can be determined. We modified the 
formulation to include the phantom field case. The equation can be simplified to linear type if we consider the fiat 
universe case fc = or the cases n = 2orn = 4 However, even if the equation is linear, it can not be considered 
as an analog to non-relativistic time-independent quantum mechanics because in this work, the wave function of 
Schrodinger-type formulation is found to be, in general, non-normalizable. Afterward, we consider a particular case 
of power-law expansion of scale factor. We show relations between cosmological quantities in conventional form and 
in Schrodinger-like form for power-law expansion. We obtain scalar field potential V{t), Schrodinger potential P{x) 
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FIG. 3: V{t) obtained from non-linear Schrodinger-type formulation for closed, flat and open universe in presence of dust and 
scalar field. 




FIG. 4: P{x) plotted versus x for power- law expansion. We set q = 2,K. = l,D — 1 and r = 0. The scalar field dominant case 
can not be plotted because even tough we set a condition D = 0, the coefficient n of the barotropic fluid equation of state still 
appears in the first and second terms of the Eq. (|36|) . There is only a real-value P{x) for the cases k = ±1 with n — 4 because, 
when X > 0, P{x) becomes imaginary in these cases. 



and wave function u{x). In the case of a scalar field dominant in flat universe, our analytical results V{4>) and 
agree well with the well-known results in A range of plots in various cases of closed, flat or open geometries is 
presented. Wave functions for the power-law expansion case (seen in the Fig. [5]) are found to be all non-normalizable 
as conjectured. 

Without knowledge of a(t), one might wonder if we could start the calculation procedure from solving the 
Schrodinger-type equation (|12p for example, the linear case as done in basic quantum mechanics. However, in order 
to do this, we must know the Schrodinger potential P{x) (Eq. ((H])) which depends explicitly on a{t) and <j>. Never- 
theless, 4> (Eq. (fTO|) also depends on a{t). Therefore we need to know the law of expansion a{t) before proceeding the 
calculation. Knowing a{t) enables us to know u{x) directly (see Eq. (|28|) ). Hence in Schrodinger-type formulation, 
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FIG. 5: u{x) plotted versus x for power-law expansion with q = 2. We set r = 0. The wave function is plotted for n = 
(cosmological constant), n — 2, n = 3 (dust), n = 4 (radiation) and n = 6 (stiff fluid). There is no real- value wave function for 
n = 3, n = 4 and n = 6 unless x < 0. 

we do not work as in basic quantum mechanics in which major task is to solve the Schrodinger equation for u{x). 
There could be many solutions of a Schrodinger-type equation. In quantum mechanics vahd solutions u{x) must be 
only normalizable type. Here, unlike in quantum mechanics, our u{x) must be non-normalizable. 

At late time the scalar field dark energy and cold dark matter (dust) are two major components of the universe 
while the others are negligible. For power-law expansion, the procedure is suitable for studying the system of scalar 
field dark energy and dark matter because it gives all real- value of P{x) for any k. We need to know a(t), k and 
D which are observable in order to find V(t). Information of V{4>) is important because it is a link to fundamental 
physics. If one starts from fundamental physics with a particular potential V{(j>) and if also knowing how (p evolves 
with t, then V could be expressed as function of t. Finally, the potential V{t) obtained from observation and another 
V{t) proposed by fundamental physics can be compared. The non-linear Schrodinger-type formulation might provide 
an alternative mathematical approach to problem solving in scalar field cosmology. 

Acknowledgements 

B. G. is a TRF Research Scholar under a TRF-CHE Research Career Development Grant of the Thailand Research 
Fund and the Commission on Higher Education of Thailand. This work is also supported by Naresuan Faculty of 
Science Research Scheme. 



[1] S. Masi et ai, Prog. Part. Nucl. Phys. 48, 243 (2002) [arXiv: astro-ph/0201137]; C. L. Bennett et ai, Astrophys. J. Suppl. 

148, 1 (2003) [arXiv: astro-ph/0302207]; D. N. Spergel et al. [WMAP Collaboration], Astrophys. J. Suppl. 148 (2003) 

175 [arXiv: astro-ph/0302209]. 
[2] R. Scranton et al. [SDSS Collaboration], [arXiv: astro- ph/0307335]. 

[3] A. G. Riess et al. [Supernova Search Team Collaboration], Astron. J. 116, 1009 (1998) [arXiv: astro-ph/9805201]; S. Perl- 
mutter et al. [Supernova Cosmology Project Collaboration], Astrophys. J. 517, 565 (1999) [arXiv: astro-ph/9812133]; 
A. G. Riess, arXiv: astro-ph/9908237; G. Goldhaber et al. [The Supernova Cosmology Project Collaboration], arXiv: 
astro-ph/0104382; J. L. Tonry et al. [Supernova Search Team Collaboration], Astrophys. J. 594, 1 (2003) [arXiv: astro- 
ph/0305008]. 

[4] D. Kazanas, Astrophys. J. 241 L59 (1980); A. A. Starobinsky, Phys. Lett. B 91, 99 (1980); A. H. Guth, Phys. Rev. D 23, 

347 (1981); K. Sato, Mon. Not. Roy Astro. Soc. 195, 467 (1981); A. Albrecht and P. J. Steinhardt, Phys. Rev. Lett. 48, 

1220 (1982); A. D. Linde, Phys. Lett. B 108, 389 (1982). 
[5] T. Padmanabhan, Curr. Sci. 88, 1057 (2005) [arXiv: astro-ph/0411044]; E. J. Copeland, M. Sami and S. Tsujikawa, Int. 

J. Mod. Phys. D 15, 1753 (2006) [arXiv: hep-th/0603057]; T. Padmanabhan, AIP Conf. Proc. 861, 179 (2006) [arXiv: 

astro-ph/0603114]. 

[6] R. M. Hawkins and J. E. Lidsey, Phys. Rev. D 66, 023523 (2002) [arXiv: astro-ph/0112139]; F. L. Williams and 
P. G. Kevrekidis, Class. Quant. Grav. 20, L177 (2003); J. E. Lidsey, Class. Quant. Grav. 21, 777 (2004) [arXiv: gr- 
qc/0307037]; F. L. Williams, P. G. Kevrekidis, T. Christodoulakis, C. Helias, G. O. Papadopoulos and T. Grammenos, 



9 



Trends in Gen. Rel. and Quant. CosmoL, Nova Science Pub. 37-48 (2006) [arXiv: gr-qc/0408056]; F. L. Williams, Int. J. 

Mod. Phys. A 20 (2005) 2481; A. Kamcnshchik, M. Luzzi and G. Venturi, arXiv: math-ph/0506017. 
[7] J. D'Ambroisc and F. L. Williams, Int. J. Pure Appl. Maths. 34, 117 (2007) [arXiv: hep-th/0609125]. 
[8] J. D'Ambroise, Int. J. Pure Appl. Maths. 3, 417 (2008), arXiv:0711.3916 [hep-th]. 

[9] A. Melchiorri, L. Mersini-Houghton, C. J. Odman and M. Trodden, Phys. Rev. D 68, 043509 (2003) [arXivrastro- 
ph/0211522]; P. S. Corasaniti, M. Kunz, D. Parkinson, E. J. Copeland and B. A. Bassett, Phys. Rev. D 70, 083006 

(2004) [arXiv:astro-ph/0406608]; U. Alam, V. Sahni, T. D. Saini and A. A. Starobinsky, Mon. Not. Roy. Astron. Soc. 354, 

275 (2004) [arXiv:astro-ph/0311364]. 
[10] D. N. Spergel et ai, arXiv:astro-ph/0603449. 
[11] W. M. Wood-Vasey et al, arXiv:astro-ph/0701041. 

[12] R. R. Caldwell, Phys. Lett. B 545, 23 (2002) [arXiv:astro-ph/9908168]; G. W. Gibbons, arXiv:hep-th/0302199; S. Nojiri 
and S. D. Odintsov, Phys. Lett. B 562 (2003) 147 [arXiv:hep-th/0303117]. 

[13] R. R. Caldwell, M. Kamionkowski and N. N. Weinberg, Phys. Rev. Lett. 91, 071301 (2003) [arXiv:astro-ph/0302506] ; 
S. Nesseris and L. Perivolaropoulos, Phys. Rev. D 70, 123529 (2004) [arXiv:astro-ph/0410309]; J. g. Hao and X. z. Li, 
Phys. Rev. D 67, 107303 (2003) [arXiv:gr-qc/0302100]; X. z. Li and J. g. Hao, Phys. Rev. D 69, 107303 (2004) [arXiv:licp- 
th/0303093]; P. Singh, M. Sami and N. Dadliich, Phys. Rev. D 68, 023522 (2003) [arXiv:hep-th/0305110]; J. G. Hao 
and X. z. Li, Phys. Rev. D 70, 043529 (2004) [arXiv:astro-ph/0309746]; M. Sami and A. Toporensky, Mod. Phys. Lett. 
A 19, 1509 (2004) [arXiv:gr-qc/0312009l; S. Nojiri, S. D. Odintsov and S. Tsujikawa, Phys. Rev. D 71 063004 (2005); 
B. Gumjudpai, T. Naskar, M. Sami and S. Tsujikawa, JCAP 0506, 007 (2005) [arXiv:hep-th/0502191]; L. A. Urena-Lopez, 
JCAP 0509, 013 (2005) [arXiv:astro-ph/0507350]. 

[14] S. Nojiri, S. D. Odintsov and M. Sasaki, Phys. Rev. D 71, 123509 (2005) [arXiv:hop-th/0504052]; M. Sami, A. Toporensky, 
P. V. Tretjakov and S. Tsujikawa, Phys. Lett. B 619, 193 (2005) [arXiv:licp-th/0504154]; G. Calcagni, S. Tsujikawa and 
M. Sami, Class. Quant. Grav. 22, 3977 (2005) [arXiv:hep-th/0505193]; H. Wei and R. G. Cai, Phys. Rev. D 72, 123507 

(2005) [arXiv:astro-ph/0509328]; B. M. Leith and I. P. Neupane, arXiv:hep-th/0702002; D. Samart and B. Gumjud- 
pai, Phys. Rev. D 76, 043514 (2007) arXiv:0704.3414 [gr-qc]; T. Naskar and J. Ward, Phys. Rev. D 76, 063514 (2007) 
[arXiv:0704.3606[gr-qc]]; B. Gumjudpai, Thai J. Phys. Series 3: Proc. of the SIAM Phys. Cong. 2007 [arXiv:0706.3467 
[gr-qc]]. 

[15] F. Lucchin and S. Matarrese, Phys. Rev. D 32, 1316 (1985). 

[16] Z. H. Zhu, M. Hu, J. S. Alcaniz and Y. X. Liu, Astron. and Astrophys. 483, 15 (2008) 
[17] C. Rubano and J. D. Barrow, Phys. Rev. D 64, 127301 (2001). 



